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The distribution of linearly polarized gluons inside a large nucleus is studied in the framework of the 
color glass condensate. We find that the Weizsacker- Williams distribution saturates the positivity 
bound at large transverse momenta and is suppressed at small transverse momenta, whereas the 
dipole distribution saturates the bound for any value of the transverse momentum. We also discuss 
processes in which both distributions of linearly polarized gluons can be probed. 

1 Introduction 
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Recently, transverse momentum dependent parton distributions (TMDs) [lj inside a nucleon have 
attracted a lot of interest. Since TMDs not only depend on the longitudinal momentum fraction x of 
the parton but also on its transverse momentum k±, they contain more detailed information on the 
internal structure of the nucleon as compared to the conventional collinear parton distributions. So 
far, the main focus of the held has been on quark TMDs. In particular, the (naive) time-reversal odd 
quark Sivers function [2] and Boer- Mulders function [3j, which are intimately linked to initial/final 
state interactions of the active quark [HE], have been under intensive investigation. In comparison, 
the available studies of (polarized) gluon TMDs [6l-H2j are still rather sparse. Among them the 
distribution of linearly polarized gluons inside an unpolarized nucleon (h^ 9 in the notation of Ref. [9]) 
is of particular interest. It is the only polarization dependent gluon TMD for an unpolarized nucleon, 
and therefore may be considered as the counterpart of the quark Boer- Mulders function. However, in 
contrast to the latter, h^ 9 is time-reversal even implying that initial/hnal state interactions are not 
needed for its existence. It has been shown that this distribution, in principle, can be accessed through 
measuring, e.g., azimuthal cos 2(f) asymmetries in processes such as jet or heavy quark pair production 
in electron-nucleon scattering as well as nucleon-nucleon scattering, and photon pair production in 
hadronic collisions |10H12j . Such measurements should be feasible at RHIC, the LHC, and a potential 
future Electron Ion Collider (EIC) [I31H3]. 

It has long been recognized that the k± dependent unpolarized gluon distribution ff (also fre- 
quently referred to as the unintegrated gluon distribution) plays a central role in small x saturation 
phenomena. Due to the presence of a semi-hard scale (the so-called saturation scale) , generated dynam- 
ically in high energy scattering, f 9 (x, k±) at small x can be computed using an effective theory which is 
also known as the color glass condensate (CGC) framework (see [15H19] and references therein). There 
are two widely used k± dependent unpolarized gluon distributions with different gauge link structures: 
(1) the Weizsacker- Williams (WW) distribution [15j ; 20, 21], and (2) the so-called dipole distribution 
which appears, for instance, in the description of inclusive particle production in pA collisions [22,23j 
(see, e.g., Refs. P2JH9] for an overview). The WW distribution describes the gluon number density 
and as such has a probability interpretation, whereas the dipole distribution is defined as the Fourier 
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transform of the color dipole cross section. Very recent work has demonstrated that both types of 
k± dependent gluon distributions can be directly probed through two-particle correlations in various 
high energy scattering reactions [244125]. These studies make use of an effective TMD factorization at 
small x in the correlation limit, where the transverse momentum imbalance of, e.g., two outgoing jets 
is much smaller than the individual transverse jet momenta. Such a factorization is suggested by the 
CGC approach. 

In this Letter, we extend the calculation of f±(x, k±) to the case of h l 9 (x, k±). To be more specific, 
we compute both the WW distribution and the dipole distribution of linearly polarized gluons in the 
CGC framework. It is shown that the WW distribution saturates the positivity bound |6j at high 
transverse momenta, and is suppressed at low transverse momenta. The dipole distribution saturates 
the bound for any value of k±. Following the procedure outlined in [23J[25] we further argue that 
the WW distribution and the dipole distribution can be accessed by measuring a cos 20 asymmetry 
for dijet production in lepton nucleus scattering and for production of a virtual photon plus a jet in 
nucleon nucleus scattering, respectively. In some sense this also extends related studies [TUHTZ] to the 
small x region. 

2 Weizsacker- Williams distribution 

We start the derivation by introducing the operator definition of the WW gluon distribution inside a 
large nucleus (6j[9], 

/' 1 1 • A 

= y^/iVwl 1 '^) + ( 2^L*a ~ 1 6 ±) xh lww(. x i k ±) , (1) 

where k\ = kj_/k± (k± = \k±\). Color gauge invariance is ensured by two (future-pointing) gauge 
links in the adjoint representation. We use 

Lc = ^ e H S /-"«-A+(C-^)^ e ^/-^(CL^==co-) | (2) 

where = A%t a with (i a )fec = —if aba and f a b c denoting the structure constants of the SU(3) 
group. When the gauge links become unity by choosing the light-cone gauge with advanced boundary 
condition, the above two distributions have a number density interpretation. Note that our convention 
for h\ 9 differs from the previous literature [6l|9] by a factor k^_/M 2 , where M is the target mass. We 
also would like to mention that the violation of translational invariance inside a large nucleus prevents 
us from shifting the coordinate y to zero, which would lead to the conventional form of k± dependent 
gluon distributions. However, as shown in Ref. [25], the gluon distributions defined in Eq. ([I]) are the 
ones which can be directly related to physical observables. 

We perform the calculation of the WW gluon distributions in the CGC framework in the light- 
cone gauge by following the standard procedure (see, e.g., Ref. [T7] for an overview). By solving the 
classical Yang-Mills equation of motion with a random color source, which has only a plus component, 
the gluon field strength tensor reads 

F+\y ± ) = d + A\ Vl _) = -U(y ± )d i L aUHy±) • (3) 

All nonlinear effects are encoded in the Wilson line U^(y±) = V exp{ig d£~a(£ _ , V±)}- The quan- 
tity a satisfies the equation — V\a a {y±) = p a (y±) with p a being the color source in covariant gauge. 
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We proceed by inserting this expression into the matrix element in ([I]) and by contracting the field 
operators in all possible ways. Due to rotational symmetry and the ordering of the Wilson lines in 
the minus direction, the only allowed contraction is 



F«(t + y)F*(y)) A = ([^la^ + y)[Ul c d{a c ](y) / 

(dla^ + y)d{ a S))JukZ + y)u C a(y)) 

6(C)(TrUHz + y)U(y))\-did{r A (Z ± ) A A (jT) , (4) 



where U\ c = U ca in the adjoint representation. In the last step, we use the propagator (a a (x)ab(y)) a = 
5 a b5{x~ —y~)TA{x± — y±) \a(x~) with r^(A;j_) = l/kj_, where comes from the correlation of color 
sources generated by a Gaussian weig ht function W A [p] = exp { - \ J d?x pa ^Sf } [IS], 

(p a (x)pb(y)) A = $ 2 {x± - y±) S(x~ - y~) X A (x~) . (5) 

One can further evaluate the contraction of two Wilson lines by expanding them in powers of a, which 
leads to 

(Tr tft(£ + y)U(y)) A = {N 2 C - 1) exp | - g 2 N c [T A (0±) - T A (£ ± )] jH d(~ A A (C)} • (6) 

Collecting all the pieces one obtains 

/V 2 — 1 r - - r i 

M ww = -fca- / dy~^y±^±^ ik± ' u [-^±di^±)\XA(y-) 

xew{-ff*N c \r A (o ± )-r A (t ± )]J™ dC-A A (C~)}- (7) 

Integrating out y± and y~ we end up with 

where Sj_ = 7ri?^ is the transverse area of the target nucleus, pa = J!^- dy~ X A (y~), and = 
Q s A^ e //^ln ^ 2 . \ is the saturation scale. By appropriate projections one can now obtain both TMD 

H QCD 

gluon distributions in a large nucleus. For the unpolarized distribution one has 



x flww( x > k -0 = S 1 M W 



WW 



N*-l S± f l2 -ik x .^ 1 



N c 4^a s J ~ s ^ H 



d ^ ±e -«x-tx M_ e -— » (9) 



in full agreement with already existing calculations [20, 21 J. For the distribution of linearly polarized 
gluons we find 

xhfl w (x,k ± ) = (&&-2tf[)M$ w 



4 ,3 ^^fU-— I- co) 
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To arrive at the result in (jlOp we made use of the Bessel functions K v (x) = ^ Jl n d9 exp{ix cos 9 + 



iv9}, and the recursion relation Kv J^ = — ^ v ^ x> ■ Note that both ff W w and h^ g ww depend, in 
particular, also on the CGC parameter Q s . 

Let us now discuss the expression in Eq. (]10p in the limit of high and low transverse momenta. For 
k± ^> Q s , the integral is dominated by small distances £j_ <C 1/Q S and can be evaluated by expanding 

the exponential exp < ^-^ >, leading to 

xhH w (x,k ± ) ~ 2S ± ^-^^ (fcj_»Q.). (11) 

For Aqcd ^ k± <^ Q s the dominant contribution comes from large distances £j_ 3> l/Q s , where one 
can neglect the exponential. We further neglect the logarithmic ^-dependence in the saturation scale 
Q s and arrive at 

xh^ ww {x,h±) ~ 2S L ^±^ (A Q cd « k± « Q 5 ) . (12) 
On the other hand, in these limits the unpolarized gluon distribution takes the form 

xfl ww (x,k ± ) ~ S± ^-^^ (fc]_»Q.), (13) 

iV 2 - 1 1 Q 2 
- fti ^ 3 ~ Ingr (Aqcd « fc± « Q s ) ■ (14) 

From those results one immediately finds that for large k± the distribution of linearly polarized gluons 
saturates the positivity limit, which in our notation reads h^ 9 < 2/f [6]. This is actually not a 
very surprising result because, like for the unpolarized gluon distribution, the correct perturbative 
tail |27H29| can be recovered for h± 9 at large k±, for which one also finds complete linear polarization. 
In contrast, the ratio h 1 9 vw / ffww is suppressed in the region of small k±, where gluon re-scattering 
effects play a more important role. 



K v + l( X ) M n t„ +V,o+ Krvt-^ f9 onrl h ± 9 



3 Dipole distribution 

We now proceed to the calculation of the dipole distribution. In that case the operator definition 
reads 



= y x/i V(^^) + (2^^"4^t^ p(x ' A;±) ' (15) 

where U [ f ] = U n (0, -co; 0)£/ n (-oo, and U [ ^ ] = C/ n (0, +co; 0)?7 n (+co, £"; £ ± ) are gauge links 

in the fundamental representation. In covariant gauge, the only nontrivial component of the field 
strength tensor is F +l (y±) = —d]_a(y±), which can be viewed as the realization of the eikonal ap- 
proximation in the McLerran-Venugopalan model. By noticing this fact, one may easily see that 
U n [+oo,^-;U]F +i (OU n H-oo,r;U] oc dyU^i-oo, +oo; [25]. This relation finally leads to 



• _ k^kjN r^± iE± . € _^L ( 

Mdp ~ 2^a s b± J (2iry & 
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where Q 2 = a s Cjr /x,4 In ^ .\ is the quark saturation momentum. Contracting M l A p with the 

q ^± QCD 

different tensors one readily finds 

xh^ Dp (x, k ± ) = 2xfl Dp (x, k ± ) = J |^ e~^<- , (17) 

which means that the positivity bound is saturated for any value of k±. Note that both ff^p and 

^idp depend, in particular, also on the CGC parameter Q sq . At large k±, the correct perturbative 
tails are recovered for both the WW and the DP distributions in the unpolarized and the polarized 
case, while the DP-type distributions are more suppressed than the WW-type distributions at small 
k±. For more discussion and additional physical insights about the difference between the two type 
distributions, see pjJJ[25] and references therein. 

4 Observables 

The extraction of gluon TMDs through two-particle correlations in various high energy scattering 
processes at small x relies on an effective TMD factorization valid in the correlation limit [241125] . 
where the transverse momentum imbalance between two final state particles (or jets) is much smaller 
than the individual transverse momenta. Normally, higher twist contributions at small x are equally 
important as the leading twist contribution because of the high gluon density. Therefore, in order to 
arrive at the mentioned effective TMD factorization, an analysis including all higher twist contributions 
would be crucial. For the unpolarized case it has been shown that the results from the effective TMD 
factorization are in agreement with the results obtained by extrapolating the CGC calculation to the 
correlation limit |24|,l25j . By applying a corresponding power counting in the correlation limit, we 
find a complete matching between the effective TMD factorization and the CGC calculation in the 
polarized case as well. For simplicity, we will express our results in terms of k±_ dependent gluon 
distributions rather than multi-point correlation functions. 

First, we discuss dijet production in lepton nucleus scattering. In fact, we consider the process 
7* + A — > q(p\) + q(p2) + X for both transversely and longitudinally polarized photons. (We also 
keep the quark mass m q in the calculation.) Because there are only final state interactions between 
the qq pair and the target nucleus, the correct gluon TMD entering the factorization formula is the 
WW distribution, in which the final state interactions, to all orders, are resummed in future-pointing 
Wilson lines. The calculation provides 

5(x^* —l)H^g^ q ql^xff ww (x,k±) 

[z 2 + (l-z q ) 2 ]e 2 P 2 -m 2 P 2 . I 

- [ Z 2 + {l-z q ) 2 ]{e) + Pi) + 2 m lPl <x»W>>Kbwl?> *±) J • (18) 

5(x r - l)H 7 *g^ q Jxfl ww (x, k ± ) + i cos{2(f))xh^ w (x, j , (19) 

where x 7 * = z q + z q , with z q , z q being the momentum fractions of the virtual photon carried by 
the quark and antiquark, respectively. The phase space factor is defined as dP.S. = dyidy2d 2 Pj_d 2 kj_, 
where yi, y<i are rapidities of the two outgoing quarks in the lab frame. Moreover, P± = (pi± —p2±)/2, 
and e 2 = z q (l — z q )Q 2 + m 2 . The transverse momenta are defined in the 7* A cm frame. In the 
correlation limit, one has |Pj_| ~ \pi±\ — \P2±\ ^> \k±\ = \pix +P2_l|- The (azimuthal) angle between 
and P± is denoted by 4>. The hard partonic cross sections i? 7 * L g^, qq can be found in Ref. [25]. Our 
calculation for these coefficients agrees with the results of [25] ■ The cos(2(/>)-modulation of the cross 



d(jl * T A^qq+X 

dP.S. 
dP.S. 
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section allows one to address the distribution of linearly polarized gluons. For intermediate values 
of x this was already pointed out in Ref. [11]. Our calculation indicates that the largest azimuthal 
asymmetry can be expected for longitudinally polarized photons. 

Let us now turn to the dipole distribution at small x. From a theoretical point of view, the simplest 
process to address h± 9 DP seems to be back-to-back virtual photon plus jet production in pA collisions, 
i.e., p + A — >■ j* (p\) + q{pi) + X. For this reaction, the multiple gluon attachments to the initial and 
final state quark line can be resummed which leads to the dipole type gluon distribution. Moreover, 
in the forward rapidity region of the proton, one may simplify the calculation by adopting an hybrid 
strategy [23J in which the dense target nucleus is treated as color glass condensate, while on the side 
of the dilute projectile proton one uses ordinary integrated parton distributions. Even though, to the 
best of our knowledge, a general proof of this method is still missing we use it here for the process 
under discussion. The differential cross section, obtained in the effective TMD factorization, reads 

dps = __ X pfl( X p) {#?9->7*<Z X fl,Dp( X > + cos ( 2 ^) H ?g S Srq xh i 9 Dp( x i k -L) 

= __ x pfi( x p) x fi,Dp( x i k ± )H gg ^ r Jl + cos(20) ^ ^ ' - j , (20) 

where the partonic cross sections are given by 

a s a em e 2 q f § u 2Q 2 i\ „ C os(2<£) _ OL s a em e 2 q / -Q 2 i 



«_->7*_- N§2 {us §&)' H ^ri- M S 2 \ - J' 
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Here we used the partonic Mandelstam variables s = (p\ + f_) 2 , u = (pi — p) 2 and t = (p2 — p) 2 , 
with p denoting the momentum carried by the incoming quark from the proton. Note that the cos 2(f) 
modulation drops out for prompt (real) photon production. In the second line in Eq. (I20p we made 
use of the relation (|17p . implying that the azimuthal dependence of the cross section is completely 
determined by kinematical factors. 



5 Summary 

We derived both the WW distribution and the dipole distribution of linearly polarized gluons in a 
large nucleus by using the CGC formalism. The WW distribution saturates the positivity bound at 
large values of while it is power-suppressed (compared to the unpolarized distribution) in the 
region of low k±. The dipole distribution saturates the bound for any value of kj_. It is worthwhile 
to point out that in the quark target model, treated to lowest nontrivial order in perturbation theory, 
h^ 9 also saturates the bound for small x and large k± [9j. We also computed h^ 9 at large k± in 
the intermediate x region within standard collinear twist-2 factorization [29]. Taking the dominant 
contribution at small x again leads to saturation of the positivity limit. It would be interesting to 
explore how the distribution of linearly polarized gluons behaves under QCD evolution effects. We 
aim at addressing this important issue in future work. 

We further argued that the WW and the dipole gluon distribution can be probed by measuring a 
cos2</> asymmetry for dijet production in DIS, and for virtual photon-jet production in pA collisions, 
respectively. Such observables can, in principle, be measured at a future Electron Ion Collider, at 
RHIC and the LHC. Studying such effects could open a new path in spin physics. Moreover, the 
results for the asymmetries constitute parameter-free predictions of the CGC framework. Therefore, 
exploring these observables may open complementary ways to test this effective theory of small x 
physics. 
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